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THE POSITIVE MASS THEOREM FOR MULTIPLE ROTATING CHARGED 

BLACK HOLES 

MARCUS KHURI AND GILBERT WEINSTEIN 


Abstract. In this paper a lower bound for the ADM mass is given in terms of the angular momenta 
and charges of black holes present in axisymmetric initial data sets for the Einstein-Maxwell equations. 
This generalizes the mass-angular momentum-charge inequality obtained by Chrusciel and Costa to 
the case of multiple black holes. We also weaken the hypotheses used in the proof of this result for 
single black holes, and establish the associated rigidity statement. 


1. Introduction 


Based on heuristic arguments reminiscent of those used to motivate the Penrose inequality (see 
Appendix]^, one may derive the following inequality 


( 1 . 1 ) 


2 ^ q^ + + 4J2 

m > --. 


relating the ADM mass m, ADM angular momentum J^ and total charge of asymptotically flat 
axisymmetric initial data for the Einstein-Maxwell equations. This inequality implies both the mass- 
angular momentum inequality m > and the mass-charge inequality m > Ig’j; the later is often 

referred to as the positive mass theorem with charge. While the mass-charge inequality has been 
rigorously established in great generality m, without the axisymmetric assumption and for multiple 
black holes, the same is not true of the mass-angular momentum inequality or the mass-angular 


momentum-charge inequality (1.1). For these inequalities, the axisymmetric condition is necessary as 


it is related to conservation of angular momentum, without which the motivating heuristic arguments 
would no longer apply. In fact, counterexamples exist [16] without the axisymmetric hypothesis. In 
this setting, and with the addition of supplementary hypotheses to be discussed below, the mass- 
angular momentum inequality was established for a single black hole by Dain in |13j . and was later 
extended and improved upon by Schoen and Zhou [22]. The case of multiple black holes was taken 
up by Chrusciel, Li, and Weinstein |9| who proved the lower bound 


( 1 . 2 ) 


m > T{Ji ,..., Jn), 


where .F is a function of the angular momentuma Jn associated with the N black holes. It is an open 
question whether this function agrees with the predicted value where J = Yln=i 


inequality (1.1) has also been settled under certain conditions for single black holes by Chrusciel and 


Costa m, in]- It is the primary purpose of the present article to extend this result to the case of 


multiple black holes, by establishing in this setting a lower bound for the mass in the spirit of (1.2). 


An initial data set (M, g, k, E, B) for the Einstein-Maxwell equations consists of a 3-manifold M, 
Riemannian metric g, symmetric 2-tensor k representing extrinsic curvature, and vector fields E 
and B which constitute the electromagnetic field. Let /tem and Jem be the energy and momentum 
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densities of the matter fields after contributions from the Maxwell field have been removed. If charged 
matter is not present, the initial data satisfy the following set of constraints 

IOtt/Iem = R + (Tr^ — \k\g — 2{\E\‘^ + 

(1.3) SttJem = divg(k — (Tr^ k)g) + 2E x B, 

divg E = divg B = 0, 

where R is the scalar curvature of g, and {E x B)i = eijiE^B^ is the cross product with e the volume 
form of g. 

It will be assumed throughout that the data are axially symmetric. This means that the group 
of isometries of the Riemannian manifold {M,g) has a subgroup isomorphic to U(l), and that all 
quantities defining the initial data are invariant under the 1/(1) action. Thus, if rj is the Killing field 
associated with this symmetry, then 

(1.4) 2,gg = £gk = £riE = = 0, 

where Tr? denotes Lie differentiation. We will also postulate that M has at least two ends, with 
one designated end being asymptotically flat, and the others being either asymptotically flat or 
asymptotically cylindrical. Recall that a domain Mend C M is an asymptotically flat end if it 
is diffeomorphic to \ Ball, and in the coordinates given by the asymptotic diffeomorphism the 
following fall-off conditions hold 

(1.5) gij = 5ij + Ol{r 2 )^ dgij £ L (Mend); kij = Ol—\{r ), /TeMj •/eM ^ ^ (Mend)) 

(1.6) K* = 0;_i(r-^), R* = Oz_i(r-^), X > ^, 
for some / > 5|^ 

Let M be simply connected. Then it is shown in [6] (see also m for the case when cylindrical 
ends are present) that M = ^^\Yln=iPn^ that there exists a global (cylindrical) Brill coordinate 
system {p, z, (j)) on M, where the points pn representing black holes all lie on the z-axis, and in which 
the Killing field is given by r/ = In these coordinates the metric takes a simple form 

(1.7) g = + dz^) + p^e~‘^^{d4> + Apdp + Azdz)‘^, 

where pe~^(dcj) + Apdp + A^dz) is the dual 1-form to and all coefficient functions are inde¬ 

pendent of (j). Let M£^^ denote the designated asymptotically flat end associated with the limit 
r=./PT^ —)• 00 . Then in this end 

11 5 3 

(1.8) 1/= Oi_3(r“2), a = o/_4(r"2), Ap = poi_^{r~2)^ A^ = 

The remaining ends associated with the points pn will be denoted by M^^, and are associated with 
the limit —>• 0, where is the Euclidean distance to pn- As the remaining ends may be either 
asymptotically flat or asymptotically cylindrical, we list both types of asymptotics 

11 13 

(1.9) [/ = 2logrn-k Oi_4(r^), a = Oi_ 4 (r,?), Ap = poi- 3 {r^), Az = oi- 3 {r^), 

(1.10) U = log rn + Oi- 4 {l), a = 0/_4(l), Ap = poi- 3 {r^), Az = oi- 3 {r^), 
respectively. 

^The notation / = oi{r~°') asserts that lim ,—>00 dif = 0 for allj < Z, and / = Oi{r~°') asserts that < C 

for all j < 1. The assumption Z > 5 is needed for the results in and m- 
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The fall-off conditions in the designated asymptotically flat end guarantee that the ADM mass, 
ADM angular momentum, and total charges are well-defined by the following limits 

1 


( 1 . 11 ) 


/ {gii,i - 9ii,jy 


( 1 . 12 ) 


Soo Soo 


(1.13) 


Qe = 


= — f 


4?! 




where S^o indicates the limit as r —)• oo of integrals over coordinate spheres Sr, with unit outer normal 
L', and A is the vector potential for the magnetic field. Due to topological considerations some care 


must be taken to construct the vector potential, moreover its contribution to (1.12) vanishes under 


appropriate asymptotic conditions HU; thus the current definition of angular momentum typically 
agrees the with the standard ADM notion. Here qe and qb denote the total electric and magnetic 
charge respectively, and we denote the square of the total charge by q^ = q^ + ql- Note that the 


fall-off in (1.5) is not strong enough to imply that the ADM linear momentum vanishes, as is typically 


assumed in the study of mass-angular momentum type inequalities. Therefore the expression (ing, 
which represents the ADM energy, does not necessarily coincide with the standard definition of ADM 
mass as the length of the 4-momentum. Nevertheless, here, we will continue to refer to ( |1.11[ ) as the 
mass. Furthermore, note that the asymptotics (1.5) are not necessarily strong enough to guarantee 


that the angular momentum is finite, since the Killing fields grow like r. However, under the addition 


hypothesis that Jem(^) £ T^(Mend) it follows that (1.12) is finite, as may be seen from the proof of 
Lemma 2.1 in |14) . 

In the presence of an electromagnetic field, angular momentum is conserved [13], mi if 


(1.14) 

In this case 

(1.15) 


J'EyAli ~ 0 - 


N 


J — ^ Jn, 


n=l 


where Jn represents the angular momentum of the black hole pn- Moreover, it will be shown in the 


next section that the condition (1.14) gives rise to a charged twist potential v which encodes the 
angular momentum by 


(1.16) 


Jn = \{v\lr -V\l^_^), 


where In denotes the interval of the z-axis between pn and Pn+i, where po = —oo and pn+i = oo. 
Potentials x s-nd may also be obtained for the electric and magnetic fields, respectively, as a result 
of the constraints div^ E = divg H = 0. Similarly, the charges of each black hole are given by 


(1.17) 

with total charges 


Qn X|/n-l)’ 


1 


Qn = - iM/n-l), 


N 

n=l 


N 


q^ = 




n=l 


(1.18) 
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In the case of a single black hole, the mass-angular momentum-charge inequality (1.1) may be 


established in two steps The first consists of proving a lower bound for the ADM mass 

in terms of a harmonic map energy functional 


(1.19) 

where 


m > M{U,v,x,'ip), 


If/ ^2u \ 

(1.20) M{U,v,X:^) = — |V[/|2 + —|Vu + xVV’-V’Vx|' + —(|VxP + |ViA| 2) dx 

87rJ^3\ p2 JJ 

with |Vf7p = {dpU)‘^ -|- {dzU)"^ and dx denoting the Euclidean volume element; this notation will 


be used throughout the paper. The inequality (1.19) relies heavily on the assumption of a maximal 


data set Tr^ A: = 0, however proposals for treating the nonmaximal case have been recently put 
forward in la, m The second step entails showing that the data arising from the extreme Kerr- 
Newman spacetime {Ukn, Vkn, Xkn, iPkn) (see Appendix [B]) , minimize the functional among all data 
with common angular momentum and charge 


( 1 . 21 ) 


X, V’) ^ AI(?7kn, i’kn, Xkn, i/^kn)- 


Since the right-hand side of (1.21) agrees with the square root of the right-hand side of (1.1), together 


with (1.19) the desired conclusion is reached. 


It should be pointed out that the hypotheses used in and |22) are unnecessarily strong. 

In these works it is assumed that the matter density is nonnegative /Uem > 0, the current density 
vanishes IJemI^ = 0, and that the 4-currents for the electric and magnetic fields (sources for the 
Maxwell equations) vanish. The later assumption concerning the 4-currents is imposed in order to 
secure the existence of potentials for the Maxwell held, and | JemIs = 0 is used to obtain a charged 
twist potential. Note that the use of 4-currents in general requires reference to an axisymmetric 
spacetime, as opposed to the initial data alone. This is justihed, since in electrovacuum the existence 
of an axisymmetric evolution of the initial data follows from its smoothness i, 0- For our purposes, 
however, reference to the spacetime can be avoided since we will show that the potentials arise in a 
direct manner from the initial data, under the weakened hypotheses div^ E = div^ B = Jem(i/) = 0. 

Theorem 1.1. Let {M,g,k,E,B) be a smooth, simply connected, axially symmetric, maximal initial 
data set satisfying plem > 0 and JEMid) = 0, and with two ends, one designated asymptotically flat 
and the other either asymptotically flat or asymptotically cylindrical. Then 


( 1 . 22 ) 


2^9^ + + 4J2 

m > --. 


and equality holds if and only if {M, g, k, E, B) is isometric to the canonical slice of an extreme 
Kerr-Newman spacetime. 

We point out that the rigidity statement of this result does not seem to have been properly 
established in the literature, even in the uncharged case. What has been previously established [22], 
is that in the case of equality the map into complex hyperbolic space arising from the given data 
agrees with the extreme Kerr-Newman harmonic map. 


In the case of multiple black holes, the first step leading to (1.19) may be established using the 


same arguments as those in the single black hole case. Thus, it is in the second step (1.21) where the 


significant difference occurs. Here the minimizing harmonic map no longer arises from the extreme 
Kerr-Newman solution, or any other well known black hole solution in general. An exception happens 
in the special situation when all charges have the same sign and the angular momenta vanish, in 
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which case the minimizing harmonic map arises from the Majumdar-Papapetrou solution. In the 
generic case, a solution (C/q, UO) Xo, i/'o) to the harmonic map equations is constructed which has 
similar asymptotic behavior to that of the extreme Kerr-Newman map near each puncture pn and 
at the designated asymptotically flat end. This asymptotic behavior allows an application of the 
convexity arguments in |22], showing that the constructed solution minimizes the functional A4 and 
yields a gap bound (see Theorem 4.1). Let 

(1.23) ■ ■ ■ ,JN,qi, • • • ■ - iQn) = ■M{Uo,vo,xo,'^o) 

denote the minimum value of the functional. Our main result is as follows. 


Theorem 1.2. Let {M,g,k,E,B) be a smooth, simply connected, axially symmetric, maximal initial 
data set satisfying Pem > 0 and Jem{i]) = 0, and with + 1 ends, one designated asymptotically flat 
and the others either asymptotically flat or asymptotically cylindrical. Then 

(1.24) m > F{Ji, q\[). 

The functions F appearing in (1.2) and (1.24) agree when the charges vanish. Hence, Theorem 


1.2 generalizes the result of [U] by including charge, and slightly improves this previous result in that 


the asymptotic assumptions on k have been weakened. Whether or not the right-hand side of (1.24) 


agrees with the square root of the right-hand side of (1.1) is an important open question. Note that 


the case of equality is not addressed in Theorem 1.2, and is closely related to the existence question 
for multiple rotating black hole solutions to the axisymmetric stationary electrovacuum Einstein 
equations. In fact we will present arguments, based on the mass gap bound of Theorem 4.1, which 


suggest that generically equality cannot be achieved in (1.24) when N > 1. 

Conjecture 1.3. Under the hypotheses of Theorem 1.2. equality in (1.24) cannot be achieved if 


N > 1 unless all charges are of the same sign and the angular momenta vanish. In this special case, 
the initial data set is isometric to the canonical slice of a Majumdar-Papapetrou spacetime. 

This paper is organized as follows. In the next section we describe a deformation of the Maxwell 
field suited for the existence of potentials, and prove Theorem o Section will be devoted to the 
construction of a minimizer for the harmonic map functional in the case of multiple black holes, and 
appropriate estimates will be established. In Section Theorem |1.2| will be proven and arguments 
supporting Conjecture |1.3 will be given. The heuristic arguments leading to 0 will be discussed 
and extended in Appendix ^ to the case when several black holes are moving apart at high velocities. 
Lastly Appendix is included to record several formulae associated with the Kerr-Newman and 
Majumdar-Papapetrou spacetimes. 


2. Proof of Theorem 11.11 

We first describe the construction of potentials, as alluded to in the introduction. Let 

(2.1) ei = e^"“( 5 p - Apd^fl), 62 = e^"“( 5 ^ - A^d^), 63 = — 

y/9H> 

be an orthonormal frame for {M,g), with dual coframe 

(2.2) 9^ = e~^~^^dp, 9^ = e~^^°‘dz, 9^ = pe~^{dcf -P Apdp -P Azdz). 

As in [3], consider the projections {E,B) of the electric and magnetic fields to the orbit space of 
g, and let F be the associated field strength defined on the auxiliary spacetime (M x M, —dt^ + g). 
That is 

(2.3) E{ei) = E{ei), B{ei) = B{ei), f = l,2, ^(es) = ^( 63 ) = 0, 
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and 

(2.4) F{ei,dt) = Ei, F{ei,ej) = eijiB\ i,j,l = 1,2,3. 

Then 

(2.5) F(r?, •) = |r?| {B{e2)9^ - Bie,)e^) , * F{rj, •) = |r?| {Eie 2 )e^ - E{ei)e^) , 

where * denotes the Hodge star operation. It follows that 

(2.6) d{E{T], •)) = |r/|(divg 5)0^ a 0^ = 0, d{*F{r], •)) = |r?|(divg E)e'^ a 0^ = 0, 

and hence there exist potentials for the electromagnetic held such that 

(2.7) dip = F{r],-), dx = *E{r],-). 

Moreover, a calculation (Lemma 4.1 of [3]) shows that 

(2.8) d {k{r]) x rj - xdi^ + V'dx) = \v\ - X divg B + i/; div^ E) 0^ a 0^ = 0, 

yielding a charged twist potential satisfying 

(2.9) dv = k{ri) x rj — xdijj + ipdx- 

As mentioned in the introduction, the advantage of these computations is that they are made directly 
from the initial data, and do not require reference to the evolved spacetime. They also show clearly 
that the conditions divgE = divgB = Jem( 0) = 0 are necessary and sufficient for the existence of 
the desired potentials when M is simply connected. 

It should be noted that (2.5) and (2.7) imply that x V’ are constant on each interval In of the 
z-axis, and 

( 2 . 10 ) 


1 


N 


^ I EiF = y^ lira ^ [ 

4^ Js^ ^ >'"^0 47r Jq 


EiFdA 


dBr„iPn) 


EiFe^^-^rl sm9ded(p 


' dB\(pn) 


= lim — f 
47r Jgi 

1 r j' 1 

= -Yl / \v\~\^0x)e-^rnSm9d9d(j) = Y,l;ix\Ir^-x\I„-l) 

JdBdpF 

where Br{pn) denotes the ball of radius r centered at pn- Similar computations yield the expressions 
for J and in (ll.lGl) and (|1.17|). 


From (2.5) and (2.7) we also find 

„4I/-2a 

\E\l + \B\l = —^{\Vx\^ + \V^P\^). 


( 2 . 11 ) 


Furthermore from (2.9) we have 
( 2 . 12 ) 
and 


fe(ei, es) = -\v\ (e2(n) + xe2(V’) - i’e2ix)) = -\ 


(2.13) 


\-2U-a 


k{e2,e3) = \r]\ ^ (ei(n) + xei(i/>) - lAei(x)) = |r/| " (SpU + x^piA - V’^pX) 


(<9^n + xdzi’ - tpdzX) , 


It follows that 
(2.14) 




\k\] > 2 {k{ei,e3f + k{e2, e^f) = 2- 


|Vn + xVV’-V'Vx|^. 
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Recall that in Brill coordinates, the scalar curvature may be expressed simply by m,m) 


(2.15) 


2e 


-2U+2a 


R = 8AU - 4Ap,,a - 4|V[/|" - {Ap,, - A.^pY , 


where A is the Euclidean Laplacian on and Ap^z = dp+dY This leads to the following mass 
formula via an integration by parts 


(2.16) 


1 


m 


= 327 L (2- 


-2U+2a 


R + p^e-^^{Ap,z - Az,pf + 4|VC/|2) dx. 


Observe that with the help of ( |2.11[ ), (2.14), and the maximal data assumption, the scalar curvature 
may be rewritten as 

R =167rpEM + \k\^g + 2 (|E|g + |B|^) 


(2.17) 


Therefore 


=167rpEM + (l^lg - 2/c(ei, es)^ - 2 k{e 2 , e^Y) + 2 {E{eYY + B{eYY) 


„6U-2a 


AU-2a 


+ 2 - 


-^IVu + xVY - V^Vxl^ + 2^^^ (|Vx|2 + IViAP) . 


1 


m 


(2.18) 


- M{U,v,x,'ip) =:^ I (4e ^^+2" (^(es)^ + ^(es)^) + p^e '^°'{Ap^z - Az,pf) dx 

./ d 3 


+ 


1 

IGvr 


-2U+2a 


(levrpEM + {\k\‘l - 2k{ei, e^Y - 2 A:(e 2 , 63 )^)) dx. 


It should be noted that this expression holds regardless of the number of ends. In the case of two 
ends (0, HH, [22]), ( |1.21[ ) and ( |2.I8[ ) yield the inequality ( |1.22[ ) in Theorem [Hi 

Consider now the case of equality in (1.22). From (1.21) and (2.18), this implies that 


(2.19) 

( 2 . 20 ) 


A —A 

■^p,z — ■^z,p^ 


Pem — 0 , E{ez) — 71(63) — 0 , 

AI(t/, u,x,V’) =-^(t7KN,-yKN,XKN,V’KN), k{ei, Cj) = k{e 3 , 63 ) = 0 , 


{U,v,X:Y) — (CkN) Pkn, Xkn, Hn)- 
, ( (2T1 ), and (|2.19|) that all components of the Maxwell field are 


According to the gap bound in [22], a map which minimizes the functional A4 must coincide with 
the harmonic map associated with the extreme Kerr-Newman spacetime, that is 

( 2 . 21 ) 

It follows immediately from 
known and agree with those induced on the t = 0 slice of the extreme Kerr-Newman spacetime, 
{E,B) = (Ekn,71kn)- 

Now observe that from (2.17), (2.19), (2.20), and (2.21) we have R = e“^"RKN, where Rkn is the 
scalar curvature of the t = 0 slice of the extreme Kerr-Newman spacetime. Using the formula (2.15) 
produces 

(2.22) = 4AUkn - 2\VUk 

However a direct computation from extreme Kerr-Newman data yields 

(2.23) = AAUkn - 2|V^7KN|^ 


^KN| — 2Ap^2a. 


so that Ap^^a = 0. We claim that this, along with the asymptotics (1.8), imply that a = 0. Note 
that it is sufficient to show that a = 0 along the z-axis. To see this, let d G (—oo, 27r) be the cone 
angle deficiency [23] arising from the metric g at the axis of rotation, that is 


(2.24) 


2tt 


= lim 


277 • Radius 


27r — I? p-^-o Circumference p-^-o 




pe 


-u 
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Since (M, g) is smooth across the axis of rotation, the angle deficiency must vanish i? = 0, and thus 


a{0,z) = 0. Note that the integral in (2.24) is not exactly 27r • Radius, but is rather the top order 


approximation to this quantity; this is all that is needed to compute the desired limit. 

We are now in a position to show that (M, g) is isometric to the canonical slice of the extreme 


Kerr-Newman solution. By (2.19) the 1-form Apdp+A^dz is closed, and hence there exists a potential 
such that dpf = Ap and dzf = A^. Consider the change of coordinates (f) = cj)-\- f{p,z), then the 


metric (1.7) takes the form 
(2.25) 




which yields the desired result g = g'xN- Lastly, observe that (2.12), (2.13), (2.20), and q:( 0 ,z) = 0 


show that the tensor k coincides with the extrinsic curvature of the canonical extreme Kerr-Newman 
slice. This completes the proof of Theorem 0 


3. Existence of the Minimizer and Estimates 


In this section we prove the existence of a minimizer for the reduced energy (1.20), having the 
asymptotics of extreme Kerr-Newman near each of the punctures pn-, and with prescribed angular 
momenta and charges. The main tool will be Theorem 2 in [25]. We denote by T the 2 ;-axis in 
and by R the axis T minus the N punctures pn- The model map $o- \ L —which we 

construct below, is not singular Dirichlet data as defined in Definition 2 in |25|, because it does not 
satisfy condition (i). Nevertheless, this condition is not used in the proof of the theorem, and the 
only key ingredient is that the reduced energy of d>o must be finite and the pointwise tension of d>o 
must be bounded with appropriate decay at infinity, which will hold true for the map constructed 
below. Consequently, Theorem 2 in |25j can be used to conclude that there is a harmonic map 
Tq: (M^ \ T,gEuc) He which is asymptotic to he. such that distjj2 (^Oj ^o) is bounded on 
\ r, and which satisfies the desired boundary conditions on the axis. 

The complex hyperbolic space is the homogeneous Riemannian manifold where the 

metric is given by 


(3.1) ds^ = du^ + e^^{dv -|- xdip — + e^^{d)^ + di!?). 

Thus the energy density of a map ^ = (u, u, x, lA) iato He is 

(3.2) £:(!>) = |Vu|2 + e^“|Vu + yViA - lAVxl^ + e^^dVxP + |ViAP). 
Note that the metric ( |3.1[ ) is invariant under the translations 

(3.3) V V — c^jJ + bx + a, + + 


for any constants a, b, c. Furthermore, the action of these translations is transitive on any slice 
Stt = = constant}, that is, consists of a single orbit of this group action. The tension of the 

map 4> is the vector field (t“, r“, r^, r’^) on the pull-back bundle given by 

(3.4) 

= Au - 2e^“|Vu + xViA - V’VxP - 6^“ (|VxP + iViAp) , 

R = Au -|- 2Vu • Vu -|- 2(Vu -|- x^V’ “ V^Rx) ' — 2e^“(Vu -|- x^i’ — V'Rx) ' (x^X + 

= Ax + 2Vu • Vx — 2e^“(Vu -|- x^V’ “ V'Rx) ' RV'j 
= AiA + 2Vtt • ViA + 2e^“(Vu -|- x'^'4’ — V'Vx) ‘ ^X- 
Thus the tension of a map vanishes if and only if the map is harmonic. 
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Note that for each set of values Qn;9n)) ^ there is an extreme Kerr-Newman 

solution with this angular momentum, electric and magnetic charge, and there is a corresponding 
harmonic map \ T —)> In view of the transitive isometric action of the translations 

above, this map ‘^n = {un,Vn,Xn,'tpn) is only determined up to constants (an,bn,Cn), where i—)• 

Vn — Cnipn + bnXn + o,n-, Xn ^ Xn + bn, and 'i/jn + Cn, as well as a domain translation z ^ z + d. 

Thus, we can set the constant values of {vn, Xn, V’n) on the component of T' on one side of the puncture 
Pn, and the values on the other side of pn will be determined by the angular momentum, electric 
and magnetic charge of Using this freedom we obtain + 1 harmonic maps where for each 
n = 1,... the values of {vn,Xn,'ipn) agree with the values of (n, x, V’) ia the map T corresponding 
to our given initial data set on the components of T' lying to both sides of the puncture pn, and 
where the values of {vn+i,Xn+i,'4’n+i) agree with those of (n,x,'0) on the unbounded components 
of the axis T'. Furthermore, we still have an overall translation available, i.e. a translation in 
which can be applied to T. This will be used in the proof of Lemma 3.1 below. 

Let us now construct the model map <ho. For the sake of simplicity, we illustrate the construction 
when N = 2 with the help of Figure clearly the construction can be carried out with any value 
of N. First, define the map to be equal to <I >3 on B, the region outside a large ball which does 
not intersect a neighborhood of the punctures, and to be equal to on small balls Bn surrounding 
the punctures, n = 1,2. These are the dark shaded regions in Figure Next, define the map 4>o 
in narrow cylindrical tubes surrounding the components of F' joining the different BnS, the lightly 
shaded regions in FigureConsider for example 62 - We pick a smooth function 0 < A( 2 ;) < 1 which 
is 0 near Bi and 1 near B2, and define 

(3.5) $0 = (1 - A)$i + Ai2 = ((1 - A)ui + Xu2, (1 - A)ni + Xv2, (1 - A)xi + Ax2, (1 - A)V'i + AV'2) • 
Finally, we extend $0 to the remaining region U, the white region in Figure so that it is smooth. 


Lemma 3.1. The reduced energy 0/$0 = iuo,vo,xo,'4’o) is finite, the tension T(<ho) has support 
inside a bounded set, and r(^>o) is pointwise bounded. Moreover, the values of {vo,xo,'fio) agree with 
those of the given data {v,x,'f’) on each component o/F'. 


Proof. The reduced energy of the extreme Kerr-Newman harmonic map is finite, see for example [3 
mm- Thus the integral of the reduced energy density over the regions Bn is clearly finite. Also 
the integral over U is finite. Thus, it is only necessary to check the integral over Cn. For clarity we 
set n = 1. Since all the quantities we seek to estimate are geometric invariants, we can now use 
the last translation available to set all the constants = ^217-1, Xi|/i = X2I/U fi’ilh = ' 02 Ui to 
zero, where Ii is the portion of the axis between pi and p 2 . This implies that p~‘^\vn\, p~^\'>pn\ and 
P~^\Xn\, n = 1,2, are bounded in Ci. We write Uq = uq + logp and Un = Un + logp, n = 1,2. It 
follows that [/o = (1 — X)Ui + XU 2 . Thus, in Ci, we have for the reduced energy density 

(3 6) =|VUoP + ^|Vno + XoV0o - 0oVxoP + ^ (|VxoP + |V0oP) 

<|VUi|2 + IVU 2 P + C {p-^\Vvi\^ + p-^XiP + P“"|V0i|2 + • • • + |VA|2) , 

where C is a constant which depends on the supremum of \Un\, p~‘^\vn\, P~^|Xn|, P~^\'f’n\- Thus 
T'($o) is bounded over Ci, and integrating over Ci clearly gives a finite quantity. It follows that the 
integral of the reduced energy over of Cn, n = 2, 3 is also finite, and hence the reduced energy of ^0 
is finite. 

Next, consider the second claim of the lemma. Since the <!>„ are harmonic, T(<hn) = 0 on Bn, 
n = 1, 2, and r(^ 3 ) = 0 on .6. Therefore the support of t($o) is contained in U Cn. 
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Figure 1. Construction of 

For the proof of the third claim of the lemma, note that since the tension vanishes on Bn and is 
clearly bounded on it remains to check the boundedness on Cn- Once again, we focus on Ci. Since 
^n is harmonic, we have |AFn| < n = 1,2. Moreover, on Ci 

(3.7) Ano = (1 - A)AC/i + AAC/s + 2VA • V(172 - Ui) + {U 2 - C/i)AA, 
so that 

|r““| <|Auo|+2£:'(l>o) 

(3.8) <|At/i| + IAC/ 2 I + 2|VA| (£:'($!) + £'{^ 2 )) + 0|AA| + 2£\^o) 

<2(1 + |VA|) (f'($i) + £'{^ 2 )) + C|AA| + 2£'{^o), 

where C = sup(|[/i| + |172|)- It follows that |r“o| is bounded. Next observe that e^“°|V7o|, e““|xo|) 
e““|Vxo|) e“°|'0o| and e'^°\V^po\ are all bounded on Ci, since the same is true of e‘^'^"\Vvn\, e^"\Xn\, 
e“"|Vxn|, e“"|'0n| and e“"|VV’n|, n = 1,2. Using this, and e^“|Aun| < C'£i'(d>„), the proof for the vq 
component now proceeds in much the same way as for the uq component: 

g 2 no|^Co| <e2“0|Awo| +C£:'($o) 

<e2“»|Aui| + e2“»|Au2| + 2|VA| (£:'($i) + £:'(l>2)) + C|AA| + C£'{^o) 

<{C + 2|VA|) (£:'(i>i) + £'(^ 2 )) + C\AX\ + C£'i^o). 


(3.9) 
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It follows that is bounded on Ci. In a similar way, we obtain that and are 

bounded on Ci, and it follows that 

(3.10) r(l>o) = (r“°)2 + ((r»)^ + 

is bounded on Ci. 

Lastly, it is immediately apparent from the construction that the values of the potentials for the 
model map agree with those of the given data on the axis. □ 

Corollary 3.2. For any set of punctures pn on the axis L and prescribed constants Xo|/n; 

V’o|/„; n = 1,...,A^, there exists a corresponding unique harmonic map iLo = (i^Oj XOj V^o) which 
is asymptotic to ^O; CLf^d satisfies 


(3.11) l/o = uo + logp = logr„ + 0(1), uo,Xo,i/’o = 0(1), as rn0. 

Proof. The existence of Tq and the fact that it is asymptotic to <ho, follow from the main result 
in [25] combined with Lemma 3.1 In order to establish uniqueness, assume that there are two 
solutions and Ti. Since the target space has negative curvature, the distance function f{x) = 
dist]g 2 (To(x), ^i(x)) > 0 is subharmonic on \ L. From this, it follows by a maximum principle 
type argument (Proposition C.4 in |9|) that / = 0. □ 


As a consequence of the fact that 'Lq is asymptotic to the model map <I>o at spatial infinity, that 
is distjj 2 (<I>o, 'Lo) —)• 0 as r —>• oo, we obtain 

(3.12) \U(i\<C on \ U^=iS5(p„), |xo| + |xo| + liAol < C” on M^. 


It is expected [9] that more refined asymptotic fall-off estimates should hold for 'Lq, which are in line 
with (4.8 )-( |lI^ . 

In order to show that the solution 4'o minimizes the functional (1.20), we will need certain estimates 
concerning the asymptotics. 


Proposition 3.3. On \ L, the solution 4'o = (Lo, pq, Xo, lAo) satisfies the estimate 

g4i7o g2i7o (J 

(3.13) |V([/o - logp)p + —|Vuo + XoViAo - + —(|VxoP + iVi/^oP) < 

and in particular near each puncture pn it holds that 

(3.14) |Vt/o| < —, iVuo + XoVi/^o - V’oVxol < ^, |Vxo| + |Vi/)o| < —. 

P rn 

Proof. This result is analogous to (2.20) in [9|, and the proof follows in the same way. As details of 
the proof were left out in |^, we provide them here in the current and more general setting. 

Let X be a point of Euclidean distance p from the axis F, and let Bpj 2 {x) be a Euclidean ball of 
radius p/2. By the Bochner identity (equation (12) in |25|), and the fact that the target manifold 
has negative curvature, it follows that the energy density £l(^o) is subharmonic A£l(^o) A 0. 
Thus, from the mean value theorem we obtain 


f(^o)U < f 


(3.15) 


Bp/2{0 
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where = vol(i?p/ 2 ) ^ Jb denotes the average. Now observe that <5(^o) < + C /and 


3p/2 


taking the average of this inequality over Bp/ 2 {x) produces 

(3.16) f(^o)U < 

Similarly we have if'('ho) < + Cjp"^, and hence 

(3.17) f'(^o)U < 


f'(^o) + 

Bp/2(x) P 


-Bp/2 (a:) P 

It remains only to show that <f^^(^o) is of order 1/p^. 

We divide this argument into two cases: away from a neighborhood of the punctures, and in a 
neighborhood of the punctures. First assume that x is away from any puncture. Then since Uq is 
uniformly bounded, we can use the arguments of (10) and (11) in [25l pages 842-843] to obtain 


(3.18) / T'('l/o)<(4 sup |[/o|" + 2 sup \Uo\ ] I jV^p, 

'B3p/4{x) 

where is a cut-off function which is 1 on Bp/ 2 {x) and vanishes outside B^p/ 4 {x). Clearly is 

of order 1//?^ in B^pn{x), so after dividing by Yo\{Bpi 2 {x)) the desired result follows. 

Assume now that x is in a neighborhood of one of the punctures, which without loss of generality 
may be taken to be the origin. From the estimate ( |3.11 ), of Uq near the puncture, we have \Uq — 
logrj < C. Furthermore C~^r{x) < r < Cr{x) on B^p/ 4 ^{x), so that | logr — logr(x)| < logC on this 
domain. It follows that \Uo — logr(x)| < C + logC. Since logr(x) is constant on B^p/ 4 {x), we can 
replace Uq by Uq = Uq — logr(x) in the arguments above to obtain 


(3.19) 


/ ^'('I'o) < 4 sup \UQf + 2 sup \Uq\ 

'Bp/2{x) Y B3p/4{x) B3p/4{x) 


' B^p/iix) 




This completes the proof of (3.13). The estimates (3.14) follow immediately from (3.13) and (3.11). 


□ 


The previous proposition is useful because the bounds it provides in (3.13) apply globally. How¬ 


ever, more detailed estimates are possible near any compact subset of the axis F which does not 
include punctures. More precisely, iko is smooth away from the punctures and satisfies the following 
asymptotics [20] : 

(3.20) UO = Cl -kC3X0 - C2V:o + 0(/) = C4 + O(/0^), XO = C2 + O(/0^), V'o = C3+ O(/0^), 

(3.21) |V(17o - logp)| = 0(p-'+^), \AUo\ < Cp-^+^, 
for some e > 0, and 

(3.22) \\/vo\=0{p), jVxol + iVV’ol = 0(p), |Vuo-k xoVi/’o -'0oVxo| = 0(p^)- 
Lastly, we will have need of the following weighted Poincare inequalities. 

Lemma 3.4. Let /3 > ^ and h € C^. If f £ is axisymmetric and satisfies (sin“^ 0)/|0=o,-7r = 0 
then 

(3.23) 


f Hr) 


< C 


^^Hr) 

2/3-1 a AC 


Jsi<r<52 sin^^+^ 9 J&i<r<S 2 sin^^ ^ 0 

where we are using r to denote the Euclidean distance to any puncture. 
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If f ^ is an axisymmetric function with f\r=Si ,52 

(3.24) f ^ 

JSi<r<S2 ' 


= 0, and /? / 0 then 

f {drff 

Si<r<S2 


r-2/3 


Proof. The first statement slightly generalizes Proposition 2.4 in [^, but the proof there easily extends 
to this situation by keeping track of boundary terms. 

For the second statement, project everything to the {p, z) plane so that the integration takes place 
over an annulus ^((ii,(52)- Now use the fact that logr is harmonic in two dimensions, and the fact 
that f\r=Si ,52 = 0 to find 

(3.25) 0 = [ V(logr) • + 2r-^Pfdrf) . 

Ja{5i,S2) JA(Si,52} ^ ^ ^ 

From this we easily get the desired inequality. □ 


4. Convexity and Proof of Theorem 11.21 
Let ik = (m, V, X, i/^) : —)> and consider the harmonic energy on a domain C 

(4.1) Eni^)= [ iVnp + e^^lVu + xVV^-V^VxP + e^^dVxP + lViAndx. 

Jn 

If n does not intersect the rotation axis F = {/? = 0}, and we write U = u + log/o, then the reduced 
energy of the map T = (U, v, x, V’) is related to the harmonic energy of T by 

(4.2) In{'I>) = En{^)+ [ (217 - logp)^, logp, 

Jdn 

where v denotes the unit outer normal to the boundary dH and 

r p2r/ 

(4.3) Xf,('I/) = / |VC/|2 + iVu + xViA - i/'VxP + ^ (|VxP + dx. 

Jn P P 

The formula ( |4.2[ ) is obtained through an integration by parts, using the fact that log p is harmonic 
on \ F. Note that X = X^a = SttAI where Ai was introduced in Section Moreover X, which 
is referred to as the reduced energy, may be considered a regularization of E since the infinite term 
J |Vlogpp has been removed, and since the two functionals differ only by a boundary term they 
must have the same critical points. 

Let To = (mo, PO) Xo, V’o) denote the harmonic map constructed in the previous section, and let 
To = (L^O) L’O) XO) i/’o) be the associated renormalized map where Uq = uq + logp. Thus, To is a 
critical point of X. It is the purpose of this section to show that Tq realizes the global minimum for 

X. 


Theorem 4.1. Suppose that T = ([/, u,x,V’) is smooth and satisfies the asymptotics (1.8)-(1.10), 
(4.8)-(4.11) with u|r = 'PoIp; x|r = XoIp; f’lr = J’olr, then there exists a constant C > 0 such that 


(4.4) 


X(T)-X(To) >C 


distL (T, To)iia: 


1 

3 


This theorem is analogous to that of Theorem 6.1 in |22) . where the role of extreme Kerr-Newman 
is now played by the (possibly) multiple black hole solution To constructed in the previous section. 
The proof in |22] is based on convexity of the harmonic energy under geodesic deformations; such a 
property is true under general circumstances when the target space is nonpositively curved. More 
precisely, let (i, e > 0 be small parameters and set Ils,e = {S < rn ior n = 1,N] r < 2/S; p > e} 
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and As^e = B 2 /s\^s,e, where B 2 /S is the ball of radius 2/5 centered at the origin. Via a cut-and-paste 
argument, it will be shown that we may assume 'll satisfies 


(4.5) supp(f7 - Uo) C B 2 / 5 , supp(n - uq, x - Xo, i> - 1 / 0 ) C 

If lilt, t G [0,1], is a geodesic in connecting ^ and ibo (this means that for each x in the 

domain, t —)■ ^t(x) is a geodesic), then = 'I'o outside i? 2/<5 ivt,Xt,'ipt) = (^^OjXOjV’o) in a 
neighborhood of As^si so that in particular Ut = Uq + t{U — Uq) on these domains. This simple 
expression for Ut together with convexity of the harmonic energy yields 


(4.6) 


dt‘^ 


2:(^t) > 2 


|V dist][j2 ('ll, 'Iio)rdx. 


Moreover, the fact that 'I'o is a critical point implies that 


(4.7) 


4t(4.,)Uo = 0. 


Theorem 4.1 then follows by integrating (4.6) and applying a Sobolev inequality. In the remainder 

in the case of 


of this section we will justify each of these steps, following closely the strategy of 
a single black hole. Most of the effort required to establish each step consists of estimating certain 
integrals. Here, however, the techniques used for these estimates will be significantly different since 
'I'o is not known explicitly, whereas in the single black hole case 'I'o is explicit as it arises from the 
extreme Kerr-Newman solution. 

Before proceeding we record the appropriate asymptotic behavior of T. Asymptotics for U are 

and if uj = dv + x'^V' “ '4’dx then 

= p'^0{r-^) 


given in ( |1.8| ), ( |1.9[ ), ( |1.10 ) 
(4.8) 


ku 


|Vx| + iViM = pO(r 


-An 


as 


00 , 


( 4 . 9 ) |a;| = p^O{r/A^)j |Vx| + |VV’| = pO{r/A‘^) as > 0 in asymptotically flat ends, 

( 4 . 10 ) |a;| = p‘^ 0 {r/l~^), |Vx| + IVi/jI = pO(r^“^) as —)■ 0 in asymptotically cylindrical ends, 

( 4 . 11 ) |u;| = 0 (, 9 ^), |Vx| + |V?/)| = 0 (p) as p^O in ns,s- 

Note that with these asymptotics T('I') is finite precisely when A > |. Moreover, one may integrate 
along lines perpendicular to T to find 

( 4 . 12 ) Ixl, IV’I = coast +/ 9 ^ 0 (r“'’') as r —00, 

( 4 . 13 ) Ixl, liAl = const +/ 9 ^ 0 (r^“^) as —)■ 0 in asymptotically flat ends, 

( 4 . 14 ) Ixl, IV'I = const + p^O(r^“^) as —)• 0 in asymptotically cylindrical ends, 

( 4 . 15 ) Ixl, liAl = const + 0 (/ 9 ^) as p 0 in 
from which it follows that 


(4.16) 

|Vn| = pO(r-^+^) 

1 as r —)■ 00 , 

(4.17) 

Vn = pO{r/A^) as > 0 

in asymptotically flat ends, 

(4.18) 

Vn = pO{r//~^) as 0 in 

asymptotically cylindrical ends. 

(4.19) 

Vn = pO{r~^^^) as 

p —)■ 0 in 
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In order to carry out the proof of Theorem 4.1 as outlined above, we must first show that it is 


possible to approximate X('I') by replacing T with a map that satisfies (4.5). This may be achieved 
as in [23 with a three step cut and paste argument. Define smooth cut-off functions 


(4.20) 


if r < 

V>\ = <26 if I < r < |, 


0 




(4.21) 


= < 


0 if Tn < 

<1 \i 6 < rn < 26, for n = 1 
1 if r-n > 25, 


,...,iV, 


(4.22) 


= 


0 

log(p/£) 

log(vT/£ 

1 


if p < e, 
y e < p < y/e, 
ii p>y/£. 


The first step deals with the region Let 

(4.23) Fi(T) = To + ¥.1(T-To) =: {UlvlxlM). 
so that T/(T) = To on \ i? 2 / 5 . 

Lemma 4.2. lim5_^o2i(T)5^(T)) = X(T). 

Proof. Write 

(4.24) X(Fi(T))=X,<i(F/(T))+Xi<,<|(Fi(T))+X,>2(F/(T)), 

and observe that ^r< 1 (Tj (T)) —7- X(T) by the dominated convergence theorem (DCT) and since Tq 

— 5 

has finite reduced energy (F^^(T)) —)■ 0. Now write 


(4.25) Xi<,<2(F/(T)) = 


|VDi|2 + 


4<-<f P 






(|VxJP + |vv^,^|^). 


1 | 2 \ 




h 


h 


h 


We have 
(4.26) 


/ 


Ii < 2 


K?-<i 


\VU\'^ + \VUo\'^ + {U-Uof\Vpl\^ 

<4<52 

where the first two terms converge to zero by the DCT and finite reduced energy of Tq, respectively. 
For the third term we may apply Holder’s inequality and the Gagliardo-Nirenberg-Sobolev inequality 
to find 
(4.27) 


V 


(U-Uo)^\Vpl\^<[ (U-Uof 

^ <452 \ -5 5 / 




1|3 


< c 


|V(C/-[/o)|2^0. 


<»-<i 


Note that the Gagliardo-Nirenberg-Sobolev inequality applies here since U, Uq G (the Sobolev 

space of square integrable derivatives) are limits of compactly supported functions. 
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Now consider I 2 , and write 

w] =ipluj + (1 - + (v - l^o)V<y95 + (xol/’ - '!/>0 X)Vv75 

+ ¥^ 5(1 - <¥>5)[(V' - V’o)V(x - Xo) - (x - Xo)V(V’ - V’o)]- 


(4.28) 


Using ( |3.12[ ) and (|4.20|) produces 
(4.29) 


h <C 


p ^(Iwp + lwol +r ^lu-uol +r |xoV’- V’oxI 




+ IV^ - V’0p|V(x - Xo)|^ + lx - X0p|V(i/) - lAo)!^)- 

The first and second terms converge to zero by the DOT and finite red uced energy of Tq. Next, by 


(3.20) and (4.19) we have p 2 |t) — tjg] —)• 0 as p —)• 0 so that Lemma 

|V(i; - uo)^ 


to show 


3.4 


applies, together with (3.12) 


Kr< 


(4.30) 


2 r^/9^ 


<c 




j.2p2 


o4r/ 


p4(7o 


+ c 


0 , 


f fL 

|<r<| 


piUo 


(ivxi" + iviAi^) + ^(ivxor + ivi/^or) 


'f 2 p 2 


where the DCT and finite reduced energy were used in the last step. A similar argument holds for 
the fourth term on the right-hand side of (4.29), whereas the fifth and sixth terms may be directly 
estimated by terms in the reduced energy of T and Tq. It follows that I 2 —)• 0. 

Consider the first term in the integral I 3 and write 


(4.31) 
so that 

(4.32) 


VX 5 = Vxo + (x - Xo)V(^J -b ifs^ix - Xo), 




|Vxll'<c- 


i<-<l P 


nVxoP , Ix-XoP , |V(x-Xo)P 

• “T o o I 


rp2p2 


p- 


V 

The first and third terms on the right-hand side may be estimated in terms of the reduced energy. 

ins 

□ 


The same is true of the second term, after an application of Lemma 3.4 Since similar considerations 
hold for the second term in I 3 , it follows that I 3 


0 . 


Consider now small balls centered at the punctures pn- Let 

(4.33) Fsi'if) = {U, vs, xs, 'ips) 
where 

(4.34) (vs, xs, i’s) = (vo, Xo, V’o) + Ps(v - wo, X - Xo, i’ - V'o), 
so that Fs(^) = To on U^^iBs{pn)- 

Lemma 4.3. lim 5 _,.oX(F 5 (T)) =X(T). This also holds i/T = Tq outside B^js- 
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Proof. Write 

(4.35) X(F5(4/)) = ^ [Ir„<s{Fsm+l5<r„<25{Fs{^))+Ir„>2siFsm] , 


N 


and observe that by DCT 

(4.36) 

Moreover 

(4.37) 


n=l 


N 


N 


Y,^r^>2s{Fsm = Y,^r^>2s{^) ^ 


n=l 


Frr.<5{Fsm = 


n=l 


AU 


^2U 


/ \VU\^ + —\ujo\^ + —(IVxor + Ivv^oH, 


where the first term on the right-hand side converges to zero again by DCT. The second and third 
terms may be estimated by the reduced energy of (and hence also converge to zero), since the 
asymptotics (1.9), (1.10), and (3.11) imply that 

(4.38) 

near each puncture. 

Now write 

(4.39) Is<r„<25{Fs{'^)) = 


< ce^o 


AU 


AU 


' 5<rn<2S 


\vu\^+ -xNI+/ + 


'&<r„<2S P 


F<rn<2S P 


h h h 

and observe that Ii —0 by DCT. In order to estimate I 2 , write 

UJS =(PSU} + (1 - ¥?5)a;o + (p - Po)V<y95 (xoV’ - '0ox)v<y95 

ips{l - ^Ps)[{'4^ - 'i/’o)V(x - Xo) - (x - Xo)V(iA - V’o)]- 


(4.40) 


Using ([4.21|) and (|4.38|) produces 
(4.41) 


h <C 


41/ 


, Xl‘^1 + 4 

6<rn<25 \ P P 


pm piU p4U 

t- I 12 ^ I 12 ^ 

—I^^ol +^r^\v-VQ\ 

fiP r^P 


IxoiA - lAoxI' 


41/ 


+ c 


/ V'o|^|V(x - xo)|^ + lx - Xo|^|V(i/; - '(/;o)P) . 

' S<rn<2S P 


The first and second terms converge to zero by the DCT and finite reduced energy of 'I'o- 


Next, assume that pn represents a cylindrical end, so that both U,Uo ~ logr„. By (3.20) and 
3 , 


(4.19) we have p 2 |ti — -^qI —)• 0 as p —)• 0 so that Lemma 3.4 applies to yield 
r AU r .y.2 


(4.42) 


/ <C 

/5<r„<25 rip'^ 


<c 

<c 


5<rn<25 P 




/ J|V(i;-uo)|^ 

' 5<rn<2S P 


AU 


AUo 


I —i^r + ^iwop 

'5<r„<25 P P 


+ c 
^0, 


r p2U p2Uo 

/ —(|VxP + |Vi/>| 2 ) + —dVxol' + IViAol') 

S<r„<25 P P 
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where the DCT and finite reduced energy were used in the last step. A similar argument holds for 


the fourth term on the right-hand side of (4.41). The fifth and sixth terms may be directly estimated 


by terms in the reduced energy of T and Tq, since 
(4.43) 


Ix-xol ^ 


sin 9 sin 9 
To establish this, observe that by the mean value theorem 

(4.44) \{x-Xo){rn,9)\=9\de{x-Xo){rn,9')\<Crn9{\Vx{rn,9')\ + \Vxo{rn,0')\)<C9 


for some 0' < 9, where we have used (3.14) and (4.10). Performing a similar calculation based at 


9 = TT, then yields (4.43), after noting that i/), i/’o behave analogously to y, xo- Lastly, in the case that 


Pn represents an asymptotically flat end, similar computations yield the desired result. It follows 
that /2 —)• 0. 

Consider the first term in the integral I 3 and write 


(4.45) 
so that 

(4.46) 


Vx<5 = Vxo + (x - Xo)^(ps + (fs'^ix - xo), 


=21/ 


^IVx^r <C 


^S<rn<25 P J 5<rn<25 V P 


=21/0 


o2I/ 


=21/ 


|Vxo| +^-^|x-Xo| +^|V(x-Xo)|' 


rlp"^ 


P^ 


The first and third terms on the right-hand side may be estimated in terms of the reduced energy. 


The same is true of the second term, after an application of Lemma 3.4 as above. Since similar 
considerations hold for the second term in la, and it follows that /a —)• 0. □ 


Consider now cylindrical regions around the axis T and away from the punctures given by 

(4.47) = {/O < e} n {(5 < for n = 1,..., r < 2/5}, 

(4.48) y^&,e = {e < P < Ve} n {5 < for n = 1,..., A; r < 2/5}. 

Let 

(4.49) Ge^ = {U,Ve,Xe,^e) 
where 

(4.50) {Vs,Xs,l/e) = {vo,Xo,'ipo) + (/e{v - vo,X - Xo,'ip - 'ipo), 
so that ^^('I') = To on p < e. 

Lemma 4.4. Fix 5 > 0 and suppose that T = Tq on y^n=i^&{Pn), then lim£_ 5 .oT(G£(T)) = T(T). 
This also holds i/T = Tq outside B 2 /s- 

Proof. Write 

(4.51) X(G,(T)) =2c,,,(G,(T)) +Xw,,(G,(T)) +XR3\(c,,,um,e)(G's(4^)). 

Since T = Tq on u//^iBs{pn), the DCT and finite energy of Tq imply that 
(4-52) X]r3\(c^^^uvv,^^^)(G£('I')) —^X(T). 

Moreover 

(4.53) IcsJG.m = / |VC/|2 + —|u;o|2 + —(|VxoP + |Vi/oP), 

dCs,s P P 
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where the first term on the right-hand side converges to zero again by DCT. The second and third 
terms may be estimated by the reduced energy of Tq (and hence also converge to zero), since 

(4.54) \U\ + \Uo\<C on \ u^^^Bs{pn) 


by dmi ). 
Now write 


c C p^U f 2 U 

(4.55) Xw,,(G,(T))= / |VC/|2+/ —\ue\^+ —(|Vx,|2 + |VV>,|2), 

JWs,s Jm.e P P 


h h h 

and notice that /i —)• 0 by DCT. In order to estimate / 2 , write 

Ue =4’£U] -b (1 - 4>e)i^0 + {v- Vo)V(j)e + (xolp “ V’ 0 X)V(/>£ 

+ </>£(! - (t> 6 )[{ip - V'o)V(x - Xo) - (x - Xo)V(i/> - i/^o)]. 


(4.56) 


Using (4.22) and (4.54) produces 

h <C [ p"^(|a;p -b Iwol^ + (loge)"^/9"^|n - uqP + (loge)"^/0"^|xoV’ - V’oXp 


(4.57) 


'm,, 


+ IV’ - V’o|^|V(x - xo)|^ + lx - Xo|^|V ('0 - -^o)!^)- 


The first and second terms converge to zero by the DCT and finite reduced energy of Tq. The third 
term may be directly estimated with the help of (3.20) and (|4.19) 


(4.58) 


[ (loge) ‘^p ^\v-vo\^<C f (loge) ^ < C'(loge) ^0 
Jm.e Jm,e 


A similar calculation holds for the fourth term. Consider now the fifth term, and use (3.20), (4.11), 
and (|4.15|) to find 


(4.59) 


[ P V'oP|V(x-Xo)P < c / P^<C£ 

Jm.e Jm,e 


The sixth term behaves in the same way, and hence /2 —t 0. 
Consider the first term in the integral I 3 and write 

(4.60) 
so that 

(4.61) 


VXe = Vxo + (x - Xo)'^(pe + </>eV(x “ Xo), 

[ ^iVx.P <c [ p-2(|Vxo|' + (log£)- Vix - xop + |V(x - Xo)P). 


All of these terms may be estimated as above, showing that I 3 —)■ 0. □ 

By composing the three cut and paste operations defined above, we obtain the desired replacement 
for T which satishes (4.5). Namely, let 

(4.62) T,,, = G,(T 5 (F](T))). 

Proposition 4.5. Let e <C 5 <C 1 and suppose that T satisfies the hypotheses of Theorem \4-l\ Then 
satisfies (4.5) and 

(4.63) 


lim limX('I' 5 £) = T('I'). 
S —^0 s —^0 
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We are now in a position to prove the main result of this section. 


Proof of Theorem f.l, By Proposition 

§0 to § 

we have 

(4.64) 


4.5 


satisfies (4.5). Thus, if is the geodesic connecting 
To to as described at the beginning of this section, then = Uq + t{Us^e ~ Uq)- Following 






dt"^ 


dt"^' 


h 


h 


with 


(4.65) 


h = ^E^sM,e) + ^ [ {dAogp)[2{Uo + t{Us,e-Uo))-logp] 


> 2 


^<5,£ 




where convexity of the harmonic energy |22) was used in the last step, and 
l2= [ 2\V{U5,e - Uo)\^ + l6{Us,e - Uq) 


pi[Uo+t{Us,e-Uo)] 


p" 


(4.66) +/ A{Us,,-Uof -^-(iVyol' + lViAoP^ 

JA,. 


>2 


^ As,e P 

|VdisW(T5,£,To)P 




since disW (T^^^, Tq) = - Uo\ on 


It remains to show that passing ^ into the integral in (4.66) is justified. For this it is sufficient to 


dt 


show that each term on the right-hand side of the equality in (4.66) is uniformly integrable. There is 


no issue with the first term since Us^e, Uq G Consider now the second and third terms, and 

write As^i; = Bs{pn)- Uniform integrability will follow if {Us^e — Uo)‘^e°'dUs,s-Uo)^ a = 2,4 is 

uniformly bounded, since then these terms may be estimated by the reduced energy of Tq. This is 
clearly the case on as U and Uq are bounded on this region. On Bs{pn), Us^e — Uq r-u log r„ if pn 
represents an asymptotically flat end and Us^e — Uq ~ 1 if Pn represents an asymptotically cylindrical 
end. Thus, the desired conclusion follows if r“*(logr„)^ is uniformly bounded, which occurs for 


0 < to < t < 1. Since to > 0 is arbitrary, we conclude that (4.6) holds for when t G (0,1]. 


We now aim to verify (|4.7|) for Choose eo < e, <5o < <5 and write 
(4.67) 


d , d ^ 


d 


< 5 o.eo 




h 


Justification for passing ^ into the integrals, for t G (0,1], is similar to the arguments of the previous 
paragraph. Then integrating by parts, and using the Euler-Lagrange equations (B.7) satisfied by To 


together with the fact that the functionals I and E have the same critical points, produces 


(4.68) 


TV 

^3 = 0{t) - ^ / 

Ja 


2{Us,e - Uo)d,Uo - / 

n=iJdBsAPn) 


2{Us,e - Uo)dMo 
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for small t, where v is the unit outer normal pointing toward Next, using that = C/q + 

t{U5,e - Uo) and = 0 yields 


— 0{i) + 


(4.69) 


r fA[Uo+t{Us^,-Uo)\ 

/ 2VUo • V{Us,e - Uo) + A{Us,e - Uo) - J -Icuop 

P 


f p 2 [Uo+t{Us^,-Uo)] 

+ / 2(17,,, - f/o) - ■ ' (iVxoP + 

4yl5o,£o 


p 


Observe that according to the first Euler-Lagrange equation of (B.7) 


(4.70) 


/ {U 5 ,e - Uo)dMo = / 

'dBsgipn) JBsf^ipn) 


piUo 

VUo • V(t/,,, - Uo) + 2{Us,, - Uo)^\ujo\^ 




=2170 


+ / {Us,e-Uo)^{\Vxor + \Vi^on. 

^Bsgipn) P 


Note that this is justified since (3.13) implies that 


(4.71) 

/ dBr^ {Pn ) 

It follows that 


iUs,e - Uo)dpUo 


<C I p ^1 logr^l = Cr^l logr^l0 as 0. 

J dBr^ (Pn) 


d f 

hrn -im,,) = / 2Vt/o • V([/,,, - Uo) + 4(1/,,, - Uo) — \uo\ 


(4.72) 


r r p'iUo 

/ 2{Us,,-Uo)dpUo+ 2([/,,,-17o)^(|VxoP + |VV’oP). 

9CSg,eQ dCSg,eg P 


This in fact vanishes, since (4.70) holds with Bsg{pn) replaced by Csg^eo- Verification of this statement 

follows from 

(4.73) 




VUo • V(t/,,, - Uo) + {Us,e - Uo)AUo 


eo 


< / |V17,,£p + |V17oP + |A17o| —T* 0 as £q —>■ 0, 

■^Oo.eo 


which is true since 17,,, = U and Uq have finite reduced energy and |A17o| < for some e > 0 

by (3.21). Hence (4.7) holds for T,,. 

Now integrating (4.6) twice and applying the Gagliardo-Nirenberg-Sobolev inequality produces 


(4.74) 


d^{'^s,s) -2:(4'o) > 2 / |VdistH2(4'5,„4'o)| dx >C[ dist^a (4'5,£, 4'o)d2: 

"V \ / toS C 


By Proposition 4.5 lim,_^o liRL,_s.oT('h,,,) = X(4'), and thus in order to complete the proof it suffices 


to show that the limits may be passed under the integral on the right-hand side. By the triangle 
inequality, it is enough to show 


(4.75) 


lim lim / distla ('h, ,, 'l>)dx = 0. 
<5^0 £^0 /to3 “c' ’ 


As mentioned in Section the geometry of complex hyperbolic space is invariant under the 
translations v = v + bx — c'ljj, x = X + c, 'll) = iIj + b. Then using the triangle inequality and direct 
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calculation produces 
(4.76) 

distH2(4'5,£,4') 

< distH2 i{Us,e, V 5 ,e ,Xs,6^ ^S,e ), iU,VS,e,X 5 ,e^'^S,e)) + disW iiU,Vs,s, Xs,s > ^ 5 ,e ),iU,V,Xs,s^ ^S,s ) ) 


+ distjjg (([/, 

V,Xs,eP 

/ 

p,2U 

<ChU-Us,e 

H- 2' 


where vs,£ = vs,e 

+ bxs,e 

(4.77) 



M 


\U-Us,e\^<l \U-Uof. 


/R3\Bi/^ 

Since U and Uq are limits in 77^ (M^) of compactly supported functions, the Sobolev inequality implies 
that U — Uq G and hence this integral converges to zero as J —>• 0. Next, we have 


(4.78) 


pl2C/ 

|u - UAel® < 


P 


12 


N 

+ +e/ 


, 121 / 


= l-'B2s(Pn) P 


12 


\v - Uol®. 


By Lemma 3.4 and (4.54) 


A2U 


1 




12 


\v - Uol <c I - Vo 


/R3\Si/i P 


(4.79) 


<C 

<c 


I 

Jm 


\v - Vo\^ 
uo 


TT |2\ 


(|Vu|^ + |Vuo| 

R3\Bi/, P 
r p^u 4[/o 

I + " iVcuol^ 


/r3\Si/, P 


+ C 


P" 

o2U _ p21/o _ 

^(|Vx|" + IVI^P) + ^(iVxoP + iViAol^), 


7r3\Bi/, P 

since p~^\v — uo|^ is bounded. This integral then converges to zero as 5 —)• 0, as each integrand 
appears in the reduced energy. Furthermore 
r ^12C/ r 


(4.80) 


^12 


iF-uor < c 


|u — Uo|® 


■"5, 




12 


< —)• 0 as e —)• 0, 


as (3.20) and (4.19) imply that |i; — uo| < Cp^ here. 

Consider now an asymptotically cylindrical end represented by pn- By choosing constants h and c 
(used to define v) appropriately in certain domains, we may assume without loss of generality that 
X) i’j Xoi V’o vanish on the axis. Therefore we have that (3.14) implies |u —u|r| < Cr~'^p^ in B 2 s{pn)- 
Moreover (4.18) yields 

(4.81) 

and similarly 

(4.82) 


\v-vo\<C-^ on B 2 s{Pn), 


Ixl + IV'I + IXol + 


< C- 


3/2 


on B 25 {Pn 
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Next, using Lemma 3.4 produces 


021 / 


'v-vof <C 


'B2siPu) P 


12 


n I— — |6 


(4.83) 


<C 

<C 

+ c 


B25(Pu) P 

f 

B2s{Pn) 

f 

B2s{Pn) 


-^\v-vo\^ 


1 2 I— — |4 

r^‘^\v - uol 

nlO 


8l— — l4 

r°\v — uor 


(|Vu|2 + iVuoP) 




,4(7 


l^|2 


,4f/o 


pi 1^1 + 


■IcPol^ 


_ ?;nl^ 


+ c 


'B2s{Pn) P P 

r r^\v — vol"^ 
'B2siPn) P® P^ 


(IV'ol'lVxol' + IXoPlW^oP)- 


From (4.81) and (4.82) it follows that 

(4.84) r>\v-v„\^ 


I0l_ _ |4 

+ 


and hence (4.83) may be estimated by reduced energies restricted to B 2 s{pn), which converge to zero 
as (5 —)• 0. Analogous arguments hold if pn represents an asymptotically flat end. We conclude that 
(4.78) converges to zero. 


Similar computations show that the remaining integrals arising from the right-hand side of (4.76) 


also converge to zero, and therefore (4.75) holds. 


□ 


Proof of Theorem 1.2~ The asymptotic assumptions on the initial data {g,k,E,B) imply that 
{U,v,x,'if) satisfy the asymptotics (1.8)-(1.10), (4.8)-(4.11). Thus Theorem[0]applies, and Theorem 


1.2 follows from (1.19) after setting 


(4.85) 


..., TTn ^ Qn ^ ^1^' * * 5^7/) A4 (Lq, uq; Xo? i/^o) ■ 


□ 


Consider now Conjecture |1.3[ and assume that equality is achieved in ( |1.24[ ) for initial data 
{M,g,k,E,B) with N > 1 black holes. If T denotes the associated harmonic map data, then fol¬ 
lowing the proof of Theorem |1.2| yields T = Tq. Arguments in Sectionsuggest that (M, g, k, E, B) 
should then give rise to a stationary axisymmetric electrovacuum extremal black hole spacetime 
with disconnected horizon, and with g conformally flat. It is likely that this spacetime falls into 
the Israel-Wilson-Perjes class, which consists of solutions to the stationary (not necessarily axisym¬ 
metric) Einstein-Maxwell equations that are distinguished by having a conformally flat orbit space. 
Moreover, since the initial data set is maximal, it would then follow from [10] that such a spacetime 
must be the Majumdar-Papapetrou spacetime. 


Appendix A. Revisiting the Heuristic Arguments 


The heuristic physical arguments which motivate (1.1) go back to Penrose’s original derivation 
of the Penrose inequality 


Typically in such arguments, it is assumed that the end state of 
gravitational collapse is a single Kerr-Newman black hole. However, a more appropriate assumption 
for the end state is a finite number of mutually distant Kerr-Newman black holes moving apart 
with asymptotically constant velocity. This should be the result, if for instance, two distant black 
holes were initially moving away from each other sufficiently fast. We will now describe the heuristic 
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arguments for the mass-angular momentum-charge inequality in this setting. It appears that this 
has not been previously considered in the literatnre. 

Let nii, Ji, Qi denote the ADM masses, angular momenta, and total charges of the end state black 
holes. Then the total (ADM) mass, angnlar momentum, and charge of the end state is m = 

J' = ^Ji^ Q = Ylli- ^ Kerr-Newman black hole these quantities satisfy the equation jHj 


(A.l) 


2 _ A ^ + ^Jj) 

“ levr + 2 + A, 


where Ai denotes horizon area. Moreover, as a function of Ai (keeping J'i and qi fixed), the right-hand 
side is nondecreasing precisely when 


(A.2) a, > AiisJqf + AJl, 

and this inequality is always satisfied with equality only for extreme black holes. Thus, computing 
the minimum value of the right-hand side of (A.l) yields 


(A.3) 


mj > 






with equality only for extreme black holes. Let mo, Jq, qo denote the ADM mass, angular momentum, 
and total charge of an initial state. Under appropriate hypotheses, such as axisymmetry and the 
existence of a twist potential, angnlar momentum is conserved j7o = ^7 = Moreover, by 

assuming that no charged matter is present, the total charge is conserved qo = Q = Qi, and since 
gravitational waves may only carry away positive energy mo > m = ^rm. 


Lemma A.l. Let ai,bi G M and let a = Y10'i, b = Ylbi. Then 
(A.4) (a4 + 62)'/"<^(4 + 62)V4_ 

Proof. Let q = 16* | and c = then 

(A.5) \b\r^ < (E l'-■l)‘'' = (E 4) ‘'A ^ c = c. 

Hence b'^ < c^. We conclude that 

(A.6) (a- + b^) < (a‘ + c‘) ^ (af + 4)= E (“* + "b. 

□ 


Lemma A.2. Let ai, 6* G M and let a = Y^Oi, b = Y^i- Then 

(A.7) \ja^ + yja^ ^ <^\j of + \l <4 + 

Proof. By Lemma |A.1| 

(A.8) (a- + 4)‘''T(EW + 4)‘'‘)'. 

Thus, it follows that 

(A.9) \/a2 + ^ \/“? + + ^i 

□ 
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Now, let ai = qi, and bi = 2J'i, then we get 
(A.IO) 


V2m = V2'^mi > Jqf + Jqf + > y + 


Squaring both sides yields the desired result (1.1). We conclude that the heuristic arguments are 
sufficiently robust to support the mass-angular momentum-charge inequality, even for spacetimes 
with multiple black holes moving apart from one another at high velocities. 


Appendix B. The Extreme Kerr-Newman and Majumdar-Papapetrou Harmonic Maps 

First we record formulas for the extreme Kerr-Newman harmonic map. Recall that in Boyer- 
Lindquist coordinates the Kerr-Newmann metric takes the form 
(B.l) 

A —a^sin^0,.o 2asm‘^0 o . n , (r^-|-— Aa^ sin^ 0 . o „ , ,o E 


-dt -|- 


(r^ -|- — a) dtdcp + 


S S 

where 

(B.2) A = + q"^ — 2mr, 

and the electromagnetic 4-potential is given by 

Qb COS 9 


sin^ edcp"^ + -dT^ + 


S = -|- cos"' 9, 


(B.3) 


A = —(dt -|- a sm‘ 
2 -/ 


(adt -|- (r^ -|- a‘^)d(j)) , 


The event horizon is located at the larger of the two solutions to the quadratic equation A = 0, 
namely r+ = m -|- m? — a? — q^, where the angular momentum is given hy J = ma. For r > r+ it 
holds that A > 0, so that a new radial coordinate may be defined by 


(B.4) 
or rather 

(B.5) 


= ^(r - m -I- \/A), 


r = r -|- m -|- 


m? — a? — q^ 
Ar 


2 / 2 I 2 

m ^ a + q 


~ 2 2 2 
r = r-|-m, m = a +q- 


Note that the new coordinate is defined for r > 0, and a critical point for the right-hand side of (B.5) 
{m? 7 ^ + q‘^) occurs at the horizon, so that two isometric copies of the outer region are encoded on 

this interval. The coordinates (r, 9, (j)) then form a (polar) Brill coordinate system, which is related 
to the (cylindrical) Brill coordinates via the usual transformation p = rsind, 2 = rcosd. Finally, 
the harmonic map (ukn, pkn, Xkn, V’kn) : \ P —)■ E!^, Fkn = ^ikn + logp, which determines the 

extreme Kerr-Newman solution is given by 

sin^ 9{2mr — q^)' 


Ukn = - ^ log 




(B.6) 


Ukn = 'ma cos d(3 — cos^ 9) — 
qar sin^ 9 


sin"^ 9 


a{q^r — ma? sin^ 9) cos 9 sin^ 9 


Xkn — 
V'kn = 


E 

q{r^ + a?) cos9 
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(B.7) 


The Euler-Lagrange equations satisfied by this and any other harmonic map T : are given 

by 

Au - 2e^“|Vu + xViA - + |Vi/'l^) = 0, 

div [e^“(Vu + xVi/> - i/>Vx)] = 0, 
div(e2“Vx) - 2e^“Vx • (Vu + xVi/^ - i/'Vx) = 0, 
div(e2“VV’) + 2e^“VV' • (Vu + xVi/^ - i/'Vx) = 0. 

Consider now the Majumdar-Papapetrou spacetime (M x (M^ \ ds"^) with 

N 

(B.8) ds^ =+ fs, f = l + Y,—^ 

1 

n=l 

where m„ = a/ (Qu)'^ + {QnV represents the mass and total eletromagnetic charge of each black hole, 
6 is the Euclidean metric, and is the Euclidean distance to each puncture. Axisymmetry may be 
imposed by choosing the punctures pn to lie on the z-axis. Cylindrical coordinates {p, z, 4>) in 3-space 
give rise to Brill coordinates with C/mp = ~ log/, and the 4-potential is given by 


N 


(B.9) 


A = Kfdt -|- y/l — 


mn{z - Zn) 


d(p, 


0< K<1. 


n=l 


The constant k relates the electric and magnetic charges to the mass by = Km-n and = 
Vl — K^rtin- Typically the Majumdar-Papapetrou spacetime is stated without magnetic charges, 
however through a duality rotation 

(B.IO) E = {cos'd)E — B = (sin-djE-|-(cosi?).B, 

magnetic charge may be introduced so that n = cos'd. Since E = nVlogf and B = Vl — log/, 
the electromagnetic potentials are obtained from (2.7) 

(B.ll) 

so that 


dXMP = npidzfdp - dpfdz), dipup = \/l - K?p{dzfdp - dpfdz), 


(B.12) 


N 

Xmp = 

n=l 


rriniz - Zn) 


N 


V’MP = V 1 - ^ 


rriniz - Zn) 


n=l 


Lastly, since this spacetime is static there is no angular momentum, and hence Ump = 0. This, 
combined with the fact that Xmp and 'i/jmp are proportional leads to a harmonic map with a 2- 
dimensional target that is isometric to hyperbolic space, namely (ump, Ump, Xmp, iAmp) : \ T —>■ 

C where Ump = ^^mp + log p- 
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